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ABSTRACT 

The recent measurement of a 1.97 + 0.04 solar-mass pulsar places a stringent lower bound on the maximum mass of compact stars 
and therefore challenges the existence of any agents that soften the equation of state of ultra-dense matter. We investigate whether 
hyperons and/or quark matter can be accommodated in massive compact stars by constructing an equation of state based on a com- 
bination of phenomenological relativistic hyper-nuclear density functional and an effective model of quantum chromodynamics (the 
Nambu-Jona-Lasinio model). Stable configurations are obtained with M > 1.97M© featuring hyper-nuclear and quark matter in color 
superconducting state if the equation of state of nuclear matter is stiff above the saturation density, the transition to quark matter takes 
place at a few times the nuclear saturation density, and the repulsive vector interactions in quark matter are substantial. 
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1. Introduction 

The masses of neutron stars are the most sensitive among 
all their parameters to the equation of state at high densities. 
Therefore, pulsar mass measurements provide one of the key ex- 
perimental constraints on the theory of ultra-dense matter. The 
masses measured in the pulsar binaries are clustered around the 
value 1 .4 Mq and have been consider as "canonical" for a long 
time. However, in recent years mounting evidence emerged in 
favor of substantially heavier neutron stars with M < 2Mq. In 
particular, the recent discovery of a compact star with a mass of 
1.97 Mq measured through the Shapiro delay provides an ob- 
servationally "clean" lower bound on the maximum mass of a 
compact star dDemorest et alj|2010l) . 

On the theoretical side it is now well-established that the 
emergence of new degrees of freedom at high densities softens 
the equation of state of matter. For example, allowing for the hy- 
perons can reduce the maximum mass of a sequence of compact 
stars below the canonical mass of 1.4 Mq. A similar reduction 
may occur if a deconfinement to quark matter takes place, al- 
though the softening of the equation of state in this case is less 
dramatic. Thus, the observation of 2Mq mass neutron star is ev- 
idence that the ultra-dense matter in neutron stars cannot be soft, 
i.e., agents that will substantially soften the equation of state are 
potentially excluded. 

We aim to study the equation of state of ultra-dense m atter 
in the light of this recent constraint (Demo rest et alJ [2010). We 
investigate to which extent one can reconcile the non-nucleonic 
components such as hyperons and two- and three-flavor quark 
matter, along with their color superconductivity, with the exis- 
tence of neutron stars with masses 2Mq. The answer(s) to the 
question above are of fundamental importance, because if hybrid 
stars featuring quark cores surrounded by a (hyper)nuclear man- 
tle exist in nature, they could provide a unique window on the 
properties of quantum chromodynamics (QCD) at high baryon 
densities under conditions not attainable in laboratory experi- 



ments (the ultra-dense matter is in equilibrium, is charge neutral 
and in /3-equilibrium with respect to weak interactions). 

Although heavy baryons (mainly E* and A hyperons) 
were considered even before the di scovery of pulsars and their 
identi fication with the neutron stars dAmbartsumvan & Saakvanl 
1960), their emergence in the cores of neutron stars is still 
very elusiv e. Treatments b ased o n relativistic density functional 
methods dGlendennind 119851 : iGlendenn ing & Moszk owskil 
11991b [Weber 1999) predict masses that are not much larger 
than the canonical mass of a neutron star which clearly contra- 
dicts modern observations. Masses on the order of < I.SMq 
wer e obtained in non-re lativistic phen omenological mod- 
els dBalberg et al J 1 1 999t iDiapo et alJ 12010). while microscopic 
models based on hyperon-nucleon potentials, which include 
the repulsive three-bod y forces, predict l ow maximal masses 
for hy pernuclear stars (Bal do et all 1 19981 120033 IVidana et al.1 
l20ll . 

The treatment of deconfined matter at ultra-high densities 
is model-dependent. We used the Nambu-Jona-Lasinio (NJL) 
model to describe the quark matter and its color superconductiv- 
ity. The model is a non-perturbative low-energy approximation 
to QCD, which is anchored in the low-energy phenomenology 
of the hadronic spectrum. The dynamical symmetry breaking, 
by which quarks acquire mass, is incorporated in this model, 
but it lacks confinement. Our ignorance of the mechanism of 
confinement requires a free parameter in the theory, which can 
be identified with the bag constant (the latter need not be the 
same as in the MIT bag model). Furthermore, it can be elimi- 
nated in favor of a more physical quantity - the transition density 
from (hyper)nuclear to quark matter. In the low-density matter 
a candidate superco nducting phase is the two-superconducting- 
colors (2SC) phase (Ba ilin & Lovdl 19841) . which at sufficiently 
high densities transforms to the th ree-flavor color-flavor-locked 
(CFL) phase dAlford et al.l 1 1999h : our present study includes 
only these two phases, but we emphasize t hat the phase struc - 
ture of matter could be more complicated (Riis ter et al.l 12005). 
The separate problem of compact stars made of strange mat- 
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ter w ith equal number of flavors of quarks (Weber & Ne greirosl 
1201 lh will not be considered here. 

Early studies of hadron-quark phase transition and quark su- 
perconductivity within the NJL model suggested that no stable 
stars (in the CFL phase) can be obtained wit hin the standard p a- 
rameterization of this model (see, e.g., iBaldo et ail d20~0 3b)). 
Recently, NJL-model based stable sequence s of hadron-quark 
stars were obtaine d for 2SC and CFL matter dKlahn et al.i r2006: 
lAlford et alJl2007h and a three-flav or crystalline color supercon- 
ducting state dlppolito et aD l200&). The key requirement for the 
existence (and stability) of these objects was the observation 
that the nuclear equation of state must be stiff and the quark 
equation of state must be supplemented by a shift (bag con- 
stant) to enable the the quar k and nuclear equations of state 
to match (Ipp olito et alj|2008b . In the latter study stable crys- 
talline color superconducting stars with masse M ~ 2Mq were 
obtained as "twin" configurations of purely nuclear counter- 
parts. Subsequently, sequences of stars containing homogeneous 
quark matter in the 2SC and CFL phases (with maximum masses 
< 1 .8M0) were obtained in the NJL model supplemented by a 
repulsive vector inte raction ((Pagliara & Schaffn er-Bielich| [2008 ; 
Lugo nes et alJl2010h . The importance of the vector interactions 
in stiff ening the quark eq uation of state has been pointed out ear- 
lier by iKlahn et alj d2006l) . Furthermore, the emergence of twin 
configurations of color-superconducting stars was confirmed in 
complemen tary studies based on va ri ations and exten sions of the 
NLJ model (iBlaschke et alJl2010allbllAlrawaill20Toh . 

Massive hybrid configurations were also constructed within 
phenomenological parameterizations of the quark matter equa- 
tion of state motivated by the MIT bag model in combination 
with nuclear equations of state that are moderately soft; large 
masses require strong quark-quark correlations in combination 
with color superconductivity in the CFL phase ( lAlford et alj 
l2005tlWeissenborn et al.ll20Tll) . 

Below we shall adopt the point of view that the nuclear equa- 
tion of state needs to be stiff above the saturation dens ity in order 
to obtain massive hybrid stars dlppolito et al.1 12008). With this 
observation as a working hypothesis we will explore the range 
of parameters allowing for compact stars featuring hyperonic 
and/or two- and three-flavor quark matter. 



2. Models 



We start with the low-density part of the equation of state, 
where the degrees of freedom are the nucleons and hyperons. 
The nuclear equation of state, as is well known, can b e con- 
structe d starting from a number of principles, see, e.g., Weber 
dl999h and|SedrakiarJd2007). Below we will work with the rel- 
ativistic mean-field models, which are fitted to the bulk prop- 
erties of nuclear matter and hype rnuclear data to describe the 
baryonic octet and its interac tions dGlendenn ing & Moszkowski 
1 199 It lLalazissis et alJll997l) We adopt the following Walecka 



Lagrangian, which includes self-interacting cr-field 

Lb = 0"b[/(^ - gajB^fi ~ \gpBT ■ Pf,) 

B Z 

1 1 2 2 

- (m B - ga-Bcr)]ilJB + -d^crd^cr - 2 ~ m o- cr 

1 2 1 1 2 

+ 2 m ^^ - 4^" ■ Pw + ' Pf 

- jbm N (ga- N o-) 3 - ^c(go- N cr) 4 

+ ^W^-m^-iF^, (1) 

e-,fT 

where the B-sum is over the baryonic octet B 
p, n, A, X ± '°, E~'°, \pB are the corresponding Dirac fields, 
whose interactions are mediated by the cr scalar, oj^ isoscalar- 
vector and p p isovector-vector meson fields. The meson and 
the baryon masses correspond to their values in the vacuum, 
the values of nucleon-meson couplings are ga-Nlm^ — 3.967 
fm, gujN/moj = 3.244 fm, g P Nhn p = 1.157 fm for nucleons; 
the hyperon-meson couplings are obtained from these by 
multiplication by factors 0.6, 0.658, and 0.6, respectively. The 
couplings in the self-interaction terms of the cr-field are given 
by b = 0.002055 and c = -0.002651. The next-to-last term in 
Eq. ([TJ is the Dirac Lagrangian of leptons, F^ v is the energy 
and momentum tensor of the electromagnetic field; we will not 
need the explicit from of other tensors in the Lagrangian (TJJi. 
The parameters above corr espond to the NL3 parametriza- 
tion dLala zissis et aTj 1 1997b . Co mputations were made also 
with the GM3 parameterization dGlendenn ing & Moszk owskil 
119911) . and we will comment on the differences below. The 
choice of this specific parametrization was made because the 
nucleonic matter has the stiffest equation of state compatible 
with the nuclear phenomenology. The mean-field pressure of 
the (hyper)nucle ar matter can be obtained from Eq. ([TJ in the 
standard fashion dWeberlll999h . 

The high-density quark matter is described by an NJL 
Lagrangian, which is extended to include the t' Hooft interac- 
tion term (oc K) and the vector interactions (oc Gv) 

L Q = tyiyfdy, - m)ifs + Gy(iA/rV) 2 

8 

+ G s ^[(0^-A) 2 + OA''75^iA) 2 ] 

a=0 

+ G D J] [ftry s t&eakWclp] K^rys 

- K (det/tftl + y 5 )»lf\ + det y #(l - y 5 )if,]) , (2) 

where the quark spinor fields \p a a carry color a — r,g,b and flavor 
{a = it . d, s) indices, the matrix of quark current masses is given 
by m — diagf(m u ,md,m s ), A a with a = 1,...,8 are the Gell- 
Mann matrices in the color space, and Aq - (2/3)lf. The charge 
conjugate spinors are defined as if/c = Ctfr T and i^c = i// T C, 
where C = iy 2 y° is the charge conjugation matrix. The partition 
function of the system can be evaluated for the La grangia n (f2]i 
negle cting the fluctuations beyond the mean-field (Ruste ret alJ 
2005). To do so, one linearizes the interaction term keeping the 
di-quark correlations A c oc r c)%iys^ ( e a bc^ b B and quark-anti- 
quark correlations cr a oc ii a a \p a a . At zero temperature the pressure 
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is given by 



P = ^jl^J dkk 2 \ei\ +AKcr u cr d cr s 

c=l a=l 

+ X Pi-Po-B*, (3) 



where e, are the quasiparticle spectra of quarks, o>q = 
G v {QM\^ u ip u + \pl*p d \QM) and <£ = 2G v (QM\i/fliff s \QM) are 
the mean field expectation values of the vector mesons co and <p 
in quark matter, p\ is lepton pressure, po is the vacuum pres- 
sure and B* is an effective bag constant. The quark chemi- 
cal potentials are modified by the vector fields as follow p.* = 
diag f(jj. u - coo, jJd-<^o,p s - <Po)- The numerical values of the pa- 
rameters of the Lagrangian are m u d = 5.5 MeV, m s = 140.7 
MeV, A = 602.3 MeV, G S A 2 =' 1.835, KA 5 = 12.36, and 
Gd/Gs = 1. The surface tension between the (hyper)nuclear 
and quark matter is not well-known, therefore we shall adopt 
the working hypothesis that this tension is high enough to pre- 
vent the formation of mixed phases. Then, the transition from 
(hyper)nuclear matter to quark matter occurs at a certain baryo- 
chemical potential at which the pressures of these phases are 
equal. This is equivalent to the condition that pressure, P, vs. 
chemical potential, p., curves for these phases cross (Maxwell 
construction). Thus, according to the Maxwell construction of 
the deconfinement phase transition, there is a jump in the density 
at constant pressure. However, the transition density itself can- 
not be fixed, because the current NJL model does not allow us 
to fix the low-density normalization of the pressure; (this is the 
consequence of the fact that this class of models does not capture 
the confinement feature of the QCD). For this reason we intro- 
duced above an additional "bag" parameter B*, which allows us 
to vary the density at which the quark phase sets-in, thus fixing 
the density of deconfinement p tI . The transition density increases 
with B* (as well as with the vector coupling Gy). For example, 
varying B* in the range -40 MeV firT 3 to 50 MeV firT 3 we find 
variations in the transition density in the range 2.4po < Pa < 4po, 
where po is the nuclear saturation density. 

Before turning to the results of calculations, we briefly sum- 
marize the physical input of the underlying models. The masses 
and coupling in the Lagrangian ([TJ are constrained by the phe- 
nomenology of nuclear and hyper-nuclear matter; we will use 
the NL3 parameterization throughout. The parameters of the 
Lagrangian (f2| are fixed to reproduce the observables of vacuum 
QCD: the pion mass m n = 135 MeV, the kaon mass thk = 497.7 
MeV, the eta-prime mass = 957.8 MeV and the pion-decay 
constant f„ = 92.4 MeV. There remain two free parameters Gy 
and B*. The first one can be used to regulate the stiffness of the 
quark matter equation of state (Gy = gives the softest equa- 
tion of state), while B* can be varied to change the density of 
the phase transition to the quark phase (deconfinement). The 
contributions to the pressure from the condensates arises from 
two phases: the two-flavor 2SC phase, where Ai = A2 = and 
A 3 + and the three-flavor CFL phase, where A] + 0, A 2 + 0, 
and A3 + 0. The strength of the pairing field is fixed by setting 
Go = Gs, which implies stronger pairing interaction than the 
one deduced from the Fierz transformation Go — 0.75Gs ■ 



3. Results 

Figure Q] shows the equations of state of (hyper)nuclear and 
quark phases and their matching via Maxwell construction for 
several values of the parameters Gy and p tr . For low values of 
the Gy parameter there are two sequential transitions from (hy- 
per)nuclear matter to the 2SC phase and from 2SC phase to the 
CFL phase. As the onset density of the quark matter is shifted 
to higher densities, the CFL phase is expelled from the density 
range relevant for the phenomenology. For higher values of the 
Gy parameter the density jump(s) are smaller, therefore the CFL 
phase can nucleate at lower densities. This tendency eventually 
leads to a smooth transition between the phases for Gy = 0.6Gs, 
if the transition density is not too high. The t' Hooft term makes 
the transition from the 2SC to the CFL phase "smoother". 

The equations of state described above were supplemented 
by an equation of state of the low-density crust dBavm et alj 
1971). These were used as an input into the Oppenheimer- 
Volkoff equations to obtain sequences of stellar configurations. 
The resulting mass-radius relationship for massive stars is shown 
in Fig. [2] together with the largest mass meas urement to date 
M = 1.97 + O.O4M dPemorest et al.l 120101) . Masses above 
the lower bound on the maximum mass are obtained for purely 
hadronic stars; this feature is prerequisite for finding similar stars 
with quark phases. Evidently only for high values of vector cou- 
pling Gy one finds stable stars that contain (at the bifurcation 
from the hadronic sequence) the 2SC phase, which are followed 
by stars that additionally contain the CFL phase (for higher cen- 
tral densities). Thus we find that the stable branch of the se- 
quence contains stars with quark matter in the 2SC and CFL 
phases. Similar results were obtained from the GM3 parameteri- 
zation of nuclear matter if, however, the hyperons were excluded 
from the consideration. Therefore, a less stiff GM3 parameter- 
ization is qualitatively equivalent to the NL3 parametrization, 
which features the softening of the equation of state due to the 
hyperons. 

Figure[3]summarizes several results for the masses and com- 
position of compact stars as a function of the parameters of the 
model Gy and p tr . First, it shows the tracks of constant maximum 
mass compact stars within the parameter space. The decrease 
of these lines with increasing vector coupling reflects the fact 
that non-zero vector coupling stiffens the equation state. In other 
words, to obtain a given maximum mass one can admit a small 
amount of soft quark matter with vanishing vector coupling by 
choosing a high transition density; the same result is obtained 
with a low transition density, but strong vector coupling, i.e., a 
stiffer quark equation of state. For low transition densities one 
finds 2SC matter in stars, which means that weaker vector cou- 
plings slightly disfavor 2SC matter. Substantial CFL cores ap- 
pear in configurations for strong vector coupling and almost in- 
dependent of the transition density (nearly vertical dashed lines 
with 5 ~ 0.1 in Fig. |3j. Note that for a high transition density 
there is a direct transition from hyper-nuclear to the CFL phase. 
For transition densities blow 3.5po a 2SC layer emerges that sep- 
arates these phases. On the other hand, weak vector couplings 
and low transition densities produce stars with a 2SC phase only. 

In Fig. [4] we show the baryon density profiles of the 
maximum-mass stars with a transition density of p tl - = 2.5po- 
For soft quark matter Gy = we obtain quark matter only in 
the 2SC phase, which extends up to about 1/3 of the star ra- 
dius. Increasing Gy has the effect of shifting the 2SC phase to 
larger radii, whereas the CFL phase develops at the center. For 
the strongest coupling studied, all phases are present in the max- 
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Fig. 2. Mass vs radius for configurations with quark-hadron tran- 
sition density p tr = 3.4po for four values of vector coupling 
Gy/Gs = 0,0.2,0.4,0.6. The purely hadronic sequence (i.e. 
the sequence that includes nucleons and hyperons) is shown 
by black solid line. The dashed lines and the gray solid lines 
show the branches where the 2SC and CFL quark phases are 
present. The filled circles mark the maximum masses of the se- 
quences. The horizontal line show s the largest mass measure- 
ment to date (IDemorest et al. 2010). 




Fig. 3. Properties of the stars as a function the free parameters 
Gy and p a . The solid lines (blue online) show the maximum 
mass configurations realized for the pair of parameters Gy and 
Ptr. The dashed (red online) curves show the amount of CFL mat- 
ter in the configurations via the ratio S = Rcfl/R, where Rcfl 
is the radius of the CFL core, R is the star radius. The parame- 
ter space to the right from 6-0 line produces CFL stars. The 
parameter space below the dashed-dotted (green online) curve 
corresponds to stars containing 2SC matter. 



imum mass star, with the combined 2SC and CFL paired quark 
matter confined within a radius that is about the half of the 
star. Note that a similar internal structure is obtained for crys- 
talline color-superconducting stars, where quark matter is con- 
fined within a r adius of ~ 7 km for maximu m mass stars with a 
radius ~ 12 km (Knippel & Sedrakian 2009). 



4. Conclusion 

Massive neutron stars are likely to develop cores composed of 
deconfined quark matter, which should be in one of the color su- 
perconducting phases. We have constructed an equation of state 
of such matter on the basis of relativistic mean-field nuclear 
functional at low densities and effective NJL model of quark 
matter supplemented by the t' Hooft and vector interactions. 
Non-rotating spherically symmetric configuration with maxi- 
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Fig. 4. Density profiles of the maximum mass configurations 
with the quark-hadron transition density p tr = 2.5p<). The vec- 
tor coupling constant varies from zero (upper panel) to Q.6Gs 
(lower panel). The hyper-nuclear and nuclear phases are marked 
by N+HYP and N respectively. 



mum masses ~ 2Mq were studied by solving the Oppenheimer- 
Volkoff equations and were compared to the recent lower bound 
on the maximum mass of a neutron star. Our equation of state 
contains two free parameters: the transition density from (hy- 
per)nuclear matter to quark matter and the vector coupling of 
quarks. We have demonstrated that stable configurations featur- 
ing color superconducting matter can be constructed if vector 
interactions are included in the quark equation of state. A pre- 
requisite for this is a stiff nuclear equation of state. We find that 
the NL3 parameterization with hyperons or GM3 parameteriza- 
tion with nucleons only are sufficiently stiff to produce maxi- 
mum masses above the lower bound. Inclusion of quark degrees 
of freedom softens these equations of state; the required lower 
bound can still be achieved because the vector interactions in 
the quark matter can stiffen the quark equation of state. 

To conclude, relativistic hypernuclear Lagrangians, which 
predict stiff hypernuclear equations of state above saturation 
density, enable us to construct stable configurations with masses 
equal and above the measured 1.97 solar-mass star. Our con- 
figuration have "exotic" matter in their interiors in the form of 
hyperons and quark matter, of which the quark matter is color- 
superconducting in the two-flavor 2SC and/or three-flavor CFL 
phases. The choice of parameters in our models geared toward 
producing stiff equations of states above the saturation density. 
Observations of pulsars with larger masses (above 2Mq) will 
seriously challenge our current understanding of dense matter 
physics above nuclear saturation. 
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